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Abstract 

I review the constraints on the mass of gravitinos that follow from considerations 
on energy loss in stars and from Big Bang Nucleosynthesis arguments. 



1 Invited talk at the 5th Workshop on High Energy Physics Phenomenology (WHEPP-5), Pune, India, 12-26 
January 1998. 



1. STELLAR ENERGY LOSS[[L|] 



la. The Sun 



Thermal stars like our sun are systems where thermal pressure balances gravity. An 
important feature of these stars is their negative specific heat. They heat up when the system 
looses energy. Indeed, when the star gives off energy, its internal plus gravitational energy 
goes down and as a cosequence the star contracts. This is so because the gravitational energy 
itself diminishes, as required by the virial theorem. Contraction, in turn, implies a rise in 
temperature and extra nuclear burning because the internal energy has to compensate for the 
decrease of gravitational energy (virial theorem, again). Now, the added nuclear energy supply 
leads to an opposite cycle of expansion and cooling (negative specific heat) and hence to a 
lesser nuclear burning. This is a self-regulating mechanism where equilibrium is maintained 
by balancing the energy loss with nuclear fusion. 

Weakly interacting particles, if not trapped in the interior of stars, drain energy. Hence 
they lead to accelerated consumption of nuclear fuel. Models of stellar evolution tolerate 
non-standard exotic energy drain mechanisms as long as their associated luminosity does not 
exceed the solar luminosity, i.e. 



lb. Degenerate Stars 

These are systems dominated by the Fermi pressure of electrons (white dwarfs) or of 
nucleons (neutron stars). In this case the internal energy (Fermi energy) is nearly independent 
of the temperature and a loss of energy is at the expense of gravitational contraction (which 
increases the the internal Fermi energy but not the temperature). In fact the system actually 
cools. 

Here also, new hypothetical particles provide for extra energy depletion which should not 
exceed 



for white dwarfs. 

lc. Type II supernovae 

Very massive stars (much heavier than the sun) ignite all nuclear fuels up to iron. Beyond 
this point no more nuclear energy is available. When the iron core reaches its Chandrasekhar 
limit, the star enters an unstable catastrophical regime: gravitational collapse follows. 



L < 10 33 erg/s 




5 x 10 3 erg/g/s 



(2) 



1 



Gravitational collapse 

After the thermonuclear cycle is completed, Si having transformed into Fe, collapse is 
triggered by the photodissociation of Fe 56 and electron capture on nuclei. The core loses 
Fermi pressure from the electrons which no longer hold gravitational pressure and the core 
implodes. Hence, v e emission dominates during infall. As the core collapses temperature and 
density rise. For densities above 2 x 10 ll g/cm 3 neutrinos become trapped in the core. When 
the inner core reaches nuclear densities (~ 3xl0 14 g/cm 3 ) a nuclear phase transition occurs, i.e. 
bound nuclei become free nucleons. This has a dramatic effect. The core stiffens because now 
non-relativistic nucleons dominate the pressure (before, the pressure of relativistic electrons 
and electron neutrinos could not halt collapse). The inner core bounces generating a shock 
wave which is responsible for the ejection of the mantle of the star. 

The hydrodynamical core collapse happens in less than a second. The relevant dynamic 
time scales are on the order of 10 to 100 miliseconds (almost free-fall). For instance, the initial 
electron-neutrino burst lasts for about 5 ms. Finally, the proto-neutron star originated that 
way cools off to a cold neutron star in a time scale of several seconds. 

Independently of the details of collapse, to form a neutron star 

E B ~ GM 2 /R ~ 3 x 10 53 erg (3) 

have beeen released. The total luminosity in electromagnetic radiation plus kinetic energy of 
the ejecta in a supernova explosion is < 10 51 erg. Emission in gravitational radiation is at 
most 1%. THE BULK OF THE BINDING ENERGY (> 99%) IS EMITTED IN FORM OF 
NEUTRINOS. In the neutronization burst ~ 10 52 erg (10% of the total energy) are emitted. 
The rest of the energy is radiated essentially as thermal neutrinos. 

Observational data from Supernova 1987A (1MB, Kamiokande) 0] imply 

E u >2x 10 53 erg (4) 

emitted over a diffusion period of the order of ten seconds. Comparing equations (4) and (3) 
we see that any additional energy drain is allowed whenever 

L x < I0 52 erg/s (5) 

2. A CASE STUDY: BOUNDS ON GRAVITINO MASSES 

The Gravitino 

The gravitino is the spin 3/2 partner of the graviton. It is a Majorana particle with only 
transverse degrees of freedom before supersymmetry breakdown. After eating the goldstino, 
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it acquires the longitudinal degrees of freedom, i.e. the ±1/2 helicities. In some recent 
models the gravitino can be superlight. Indeed, models where gauge interactions mediate 
the breakdown of supersymmetry 0, models where an anomalous U(l) gauge symmetry 
induces SUSY breaking |3J], and no-scale models are all examples of models where a superlight 
gravitino can be accommodated ||. In all of them, the gravitino is the LSP and, furthermore, 
its couplings to matter and radiation are inversely proportional to its mass. For very small 
gravitino masses, the longitudinal components of the gravitino dominate the interactions with 
matter so that effectively, 

G, = ^|m-/ 2 9 M x (6) 

In a subclass of the models mentioned above, the scalar and pseudoscalar partners of the 
goldstino are also ultralight. Their couplings do also show the enhancing tn^j 2 factor 0. 

2a. Astrophysics M 

Supernova constraints^, |9], |10 , [TT 



According to Luty and Ponton[T2|], photons have an effective coupling to gravitinos with 
the following structure 

5L = e -{M/A 2 fdVxF,u + h.c. (7) 

where M is a model dependent supersymmetric mass parameter and A is the supersymmetry 
breaking scale. This interaction provides the leading contribution to gravitino pair emissiv- 
ity via electron-positron annihilation, nucleon-nucleon bremsstrahlung, plasmon decay, and 
photon-electron scattering followed by radiation of the gravitino pair. 

The energy-loss rate (per unit volume) via pn — > pnGG is, 

Q = [ d Lh—- ^- [j d3p * f l/2 (i- /3 )(i- /4 )(2 7 r)^ 4 (P f -P,) V | M fi | 2 {kl + k° 2 ) 
V J (2vr) 3 2A;?(27r) 3 2^/ = L 1 (2vr) 3 2p°' / U V J6JK J JK ' V ; s ~^ s 

(8) 

where {p°,p)i are the 4-momenta of the initial and final state nucleons, (k°,k)i^ are the 4- 
momenta of the gravitinos and are the Fermi-Dirac distribution functions for the initial 
proton and neutron and (1 — ^3 4) are the final Pauli blocking factors for the final proton and 
neutron. The squared matrix element can be factorised as follows, 

£ I Mfi | 2 = (2n)W(M/AYN, u G^ 2 (9) 

spins 

where is the nuclear (OPE) tensor and G^ 2 is the gravitino tensor in the matrix element 
squared. The factor is common to any bremsstrahlung process involving nucleons. It ap- 
pears, e.g., in neutrino bremsstrahlung calculations and in axion bremsstrahlung calculations, 
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and is given explicitly in ref [1|]. On the other hand, G^ 2 is a tensor specific to gravitino 
bremsstrahlung. It reads, 

G% 2 = k^k v 2 + h%k\ - k x .k 2 g^ (10) 

The integration of N^ v over the phase-space of the nucleons can be performed explicitly and 
the details can be found in Raffelt's book When we contract the result with the gravitino 
tensor G3/2 and perform the integrals over gravitino momenta to complete the energy depletion 
rate, we are led to the following emissivity: 

QZLs = (8l92/385n^)aW n (M/ArY e nlT^/mf (11) 

for non-degenerate and non-relativistic nucleons (a^ is the pionic fine-structure constant, ub 
is the number density of baryons, and Y e is the mass fraction of protons). However, nucleons 
are moderately degenerate in the SN core. The emissivity in the (extreme) degenerate case 
is calculated to be, 

Q° ems = (164tt 3 /4725)«X(M/A 2 )VT 8 (12) 

with pf, the Fermi momentum of the nucleons. Numerically, for the actual conditions of the 
star, both emissivities differ by less than an order of magnitude (about a factor of three). 
Since the actual emissivity interpolates between these two values, we shall adopt the smallest 
of the two (i.e. QZims) t° make our (conservative) estimates. We turn next to the annihilation 
process. 

The energy loss for the process e + {pi) + e~{p 2 ) — > G(ki) + G{k 2 ) can be calculated along 
similar lines as above. The spin averaged matrix element squared is in this case, 

£ I M fi | 2 = (2n) 2 a\M/A 2 ) i E, u ( Pl ,p 2 )G^ 2 (k 1 ,k 2 ) (13) 

spins 

where E^ iu (p 1 ,p 2 ) equals formally the tensor G^J 2 in eq.(10) with hi, k 2 replaced by Pi,p 2 . The 
luminosity then is found to be, 

Q ann = 8a 2 (M/A 2 ) 4 T 4 e-^ T /x 5 6(/i/T)/157r 3 (14) 

with b(y) = (5/6)e^- 5 (F 5 + F 4 " + F 4 + F 5 ") where F±(y) = / °° dxx m ~ l /(l + e x±y ) (/i is the 
chemical potential of the electrons). The function b(y) — ► 1 in the degenerate limit. Finally, 
our estimate of the plasmon decay luminosity is, 

Q P = 16C(3)a 4 T 3 /i 6 (M/A 2 ) 4 /81vr 5 (15) 

(where only transverse plasmons have been taken into account). 
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The process 7e — ► eGG has not been evaluated analitically but numerically has been seen 
to be of the same order as e + e~ — > GG^M. Taken at face value, the bremsstrahlung rate is 
the largest of the four. However, Qbrems is overestimated since we did not consider multiple 
scattering effects which are present in a dense medium[|I|. Indeed, as for the axion case[I|], 
the gravitino bremsstrahlung rate probably saturates around 20% nuclear density and this 
should be taken into account when evaluating eq.(6). If we use now the values T = 50MeV, 
H = 300MeV, and Y e = 0.3, eqs. (6) (with n B ~ 0.2n nuc ), (9) and (10) give 

Qann ' Qbrems = Qp ~ 1-2 X 10 3 : 3 X 10 2 : 1 (16) 

Therefore, a limit on A will follow from the requirement that L 3 / 2 ~ VQ ann (V is the 
volume of the stellar core) should not exceed 10 52 ergs/s. This constraint on the gravitino 
luminosity L 3 / 2 implies, in turn, 

A > 300GeV(M/A3GeV) 1/2 (T/50MeV) u/w (R c /WKm) 3/8 (17) 

or, using m 3/2 = 2.5 x \Q~ A eV(A/lTeV) 2 , 

m 3/2 > 2.3 x 10 _5 eV. (18) 

Of course, the previous calculation makes sense only if gravitinos, once produced, stream 
freely out of the star without rescattering. That they actually do so, for A > 300GeV, can 
be easily checked by considering their mean-free-path in the core. The main source of opacity 
for gravitinos is the elastic scattering off the Coulomb field of the protons: 

A = 1 1 an = (4/Tra 2 )Y- 1 p- 1 m- 1 (A 2 /M) 4 (19) 

The thermally averaged cross-section for elastic gravitino scattering on electrons is roughly 
a factor Tfi/m 2 smaller than that on protons and thus it does not contribute appreciably to 
the opacity. Putting numbers in eq.(19) we find: 

A ~ 1.4 x 10 7 cm(A3GeV/M)\A/300GeV) 8 (20) 

On the other hand, the calculation of Q breaks down for A < lOKm, i.e. for A < 220GeV, 
when gravitinos are trapped in the SN core. In this case, gravitinos diffuse out of the dense 
stellar interior and are thermally radiated from a gravitino-sphere -R 3 /2- Because in this 
instance the luminosity is proportional to T 4 , only for a sufficiently large i? 3 / 2 (where the 
temperature is correspondingly lower), the emitted power will fall again below the nominal 
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10 52 erg/s. Consequently, gravitino emission will be energetically possible, if A is small enough. 
The gravitino-sphere radius can be computed from the requirement that the optical depth 



R 



dr/X(r) (21) 
be equal to 2/3 at R = R 3 / 2 . Here, A(r) is given in eq.(19) with the density profile ansatz: 

p(r) = Pc (R c /r) m (22) 

with p c = 8 x 10 14 g/cm 3 , R c = lOKm and m = 5 — 7 and which satisfactorily parameterises 
the basic properties of SN1987A fll3"l. An explicit calculation renders: 



R 3/2 = J R c [(8F e /37ra 2 )(A 2 /M) 4 (m - \)j Pc R c m p \ x l x - m (23) 

Stefan-Boltzmann's law implies for the ratio of gravitino to neutrino luminosities, 

L 3/2 /L u = (R 3/2 /R u ) 2 [T(R 3/2 )/T(R u )} 4 (24) 

where R v is the radius of the neutrinosphere. To proceed further we use the temperature 
profile: 

T = T c (R c /r) m/3 (25) 

which is a consequence of eq.(22) and the assumption of local thermal equilibrium. Now, 



taking m = 7||14j|, we obtain 

L 3/2 /L u = {R v / R c ) 22/3 [(16Y e / ira 2 ) (A 2 /M) 4 / p c R c m p ) 11/9 (26) 

By demanding that L 3 / 2 < 0.1L U and using R v ~ 30Km, we get 

A < 70GeV. (27) 

This in turn implies m 3 / 2 < 10~ e eV. Since, on the other hand, the anomalous magnetic 
moment of the muon already requires m 3 / 2 to be larger than ~ lO -6 ^ [l^|, we are forced 
to conclude that 

A > 300GeK (28) 

or, equivalent ly, 

m 3 / 2 > 2.2 x l(r 6 eV: (29) 

In conclusion, we have carefully derived the bounds on the superlight gravitino mass (i.e. 
the SUSY scale A) that follow from SN physics. These limits are completely general in the 
sense that they do not rely on other particles in a given particular model being light. 
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Should other particles such as the scalar partners of the goldstino also be light, then the 
resulting bounds are necessarily tighter. In such clearly less general frame, constraints have 
also been derived in the literature. Before reviewing these bounds, let us refer to the very 
recent claim by Clark et al.[]I7| that the dimension 6 operator 7GG in L e ff derived by Luty 
and Ponton is actually not there. If true, the bounds derived above are 1 to 2 orders of 
magnitude worse [IT - 



Suppose now that the gravitino and S and/or P are very light (<C T). Because now these 
scalars can be emitted by astrophysical bodies, one has to consider additional energy-loss 
channels^, H|.The relevant interaction is given by 



K 

4 V 3 v m 3/2 ' 



~ l L mt = -- x l-(—L)(SF^F^ + PF, V F^) (30) 



where k = (M) l/2 Mpl . 



The main energy drain mechanism is the Primakoff process 7e — > eS/P via one pho- 
ton exchange in the t channel. Calculation of luminosities in the supernova goes along the 
same lines as before using, of course, the Primakoff cross section and which was derived in 
reference [0]. The restriction 

L < 10 52 erg (31) 



then implies 



m 



™ 3/2 > 30(— -L- )eV (32) 



3/2 - ~~ x 100GeV 
Since 

T(S/P 27) - ^(^) 2 ™!/p (33) 
967r m 3 / 2 ' 

sufficiently heavy S/P will decay inside the core. For this NOT to happen, it is required that 

m s/P < lOMeV (34) 



Also, S/P leave the supernova core without further rescattering provided 

m 3/2 > 0.3eV (35) 

(because their m.f.p. is ~ 10 10 (^^-) 2 cm). 

For masses below about 0.3 eV, S/P get trapped and their energy is radiated from a S/P- 
sphere following the Stefan-Boltzmann law L oc R^ P T^ P . The luminosity L is compatible 
again with observation (SN1987A) for m 3 / 2 < 10~ 1,5 eV. 
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Limits from the Sun^Q 

Should ms/p < IkeV then these particles could be emitted from the sun and the lower 
limits to their masses are in this case slightly different. In the sun, the m.f.p. of S/P is 

Xs/P „ io 41 (^) 2 cm (36) 



Xs/p exceeds the solar radius for m 3 / 2 > 10 3 - 5 eV 

(for rrij = lOOGeV^). The emissivity via Primakoff scattering turns out to be 

K 2 a ^fT 7 F(K D /T)V Sun (37) 
m 3/2 

where the function F(x) (see 0, p. 169) takes care of plasma effects characterized by the 
Debye momentum kd- This emissivity is bounded above by the sun luminosity, i.e. 10 33 erg/s. 
As a consequence, the gravitino mass should verify 

m 3/2 > 50eV (38) 

On the other hand, in the trapping regime (i.e. m 3 / 2 < 10~ 3 ' 5 ), S/P emission is allowed as 
long as m 3 / 2 < 10 -6 for in this case their thermal radiation is sufficiently slowed down. Since 



(g — 2)^ implies already m 3 / 2 > 10 6 |l6l we conclude 

m 3/2 > 50eV (39) 

for m 3 / 2 < IkeV and for rrij ~ O(100)GeV. 
2b. Big Bang Nucleosynthesis 

Primordial Helium-4 abundance depends on the effective d.o.f. at nucleosynthesis because 
H ~ g X J 2 Mp]T 2 and, the larger H the sooner do weak interactions decouple. Hence, the 
neutron to proton ratio is larger and the helium yield rises accordingly. 

The effective number of d.o.f. for a given species i is 
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KlZZn = Cr/sM^f (40) 



Bounds from Yp are usually presented as bounds on AN„ (i.e. number of extra equivalent 
neutrino species). It is explicitly given by the formula 

AiV,(L~) = (i /7 )|[^#\] 4/3 (41) 
1 9A L D i ) 
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where is the decoupling temperature of species i. 



The latest analysis from the group in Chicago fll9|1 sets the limit 



AN V < 1 (42) 

which implies NO bound for the superlight gravitino (equivalent to one neutrino species). 
That is, the gravitino does not have to decouple from the cosmological thermal bath prior 
to the BBN era. On the other hand, if the particles S/P are light ( ms/p < lMeV, i.e. 
relativistic at nucleosynthesis) then, they should freeze out before the Universe cools down to 
T ~ 200MeV so that AN V is very small( i.e. g m (T v ) / 'g*(T Di ) < 1). 

S/P are kept in equilibrium through 7e~ <-> S/Pe~ and this rate falls below the Hubble 
rate H when 

. 2 a(^) 2 T 3 ~ <tf>*L- (43) 
m 3/2 Mpi 

Taking T ~ 200MeV ( and = lOOGeV) we get@ 

m 3/2 > ley. (44) 

3. SUMMARY 

In a wide class of supergravity models with a supersymmetry breaking scale in the TeV 
range, the gravitino can be very light. In fact, its mass could lie anywhere between lfieV 
and IkeV. It is also a generic feature of some of the recently considered models that the 
superlight gravitino is accompanied by a superlight scalar S and pseudoscalar P particles, 
which are the neutral scalar partners of the goldstino. For energy scales such that E ^> 7713/2 
the longitudinal component of the gravitino dominates and the gravitino effectively behaves 
as a spin- 1/2 Goldstino. The neutral scalars and the gravitino are coupled to matter with 
strength inversely proportional to the gravitino mass and, hence, they can be abundantly 
produced in the interior of stellar cores. By requiring that the radiated power does not 
overcome 10 52 erg/ s in a supernova explosion and 10 33 erg/s in the case of solar emission, we 
obtain the following results. 



independently of 777-3/2- Also, 



or 



m 3/2 > !0' 5 eV (45) 
m 3/ 2 > 30eV (46) 
m 3/ 2 < OMeV (47) 
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if m s/ P < lOMeV and, 

m 3/2 > 50eV (48) 

for ms/p < IkeV. 

Finally, from Big Bang Nucleosynthesis arguments we infer the limit 

m 3 /2 > leV (49) 
should m s /p be lighter than about 1 MeV. 

/ am grateful to the organizers of the Workshop for the friendly hospitality and for the 
nice atmosphere provided in Pune. 

Work partially supported by the CICYT Research Project AEN95-0882. 
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